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Abstract 



The Lai_xSrxMn03 system with 0.2 < x < 0.4 has traditionally been 
modelled with a "double exchange" Hamiltonian, in which it is assumed 
that the only relevant physics is the tendency of carrier hopping to 
line up neighboring spins. We present a solution of the double exchange 
model, show it is incompatible with many aspects of the resistivity data, 
and propose that a strong electron-phonon interaction arising from a 
Jahn- Teller splitting of the outer Mn d-level plays a crucial role. 
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The Lai_xAxMn03 system (A represents a divalent alkali element such as Sr or 
Ca) has attracted much recent attention because of the very large magnetoresistance 
exhibited for 0.2 < x < 0.4 |1|]. Treatments of the physics of this system have 
fo cussed primarily on the phenomenon of "double exchange" 0. In this commu- 
nication we show that a straightforward and reasonably complete solution of the 
double exchange model is possible. From this solution we calculate the frequency 
and temperature dependence of the conductivity and show that it disagrees in several 
important respects with the experimental data. The discrepancy must be resolved by 
including additional physics, which we suggest involves polaron effects due to a very 
strong electron-phonon coupling coming from a Jahn- Teller splitting of the Mn'^+ ion. 
We mention several experiments which have apparently not yet been performed and 
which would verify or falsify our suggestion. 

In Lai-xAxMnOs the electronically active orbitals are the Mn dx2_y2 and d^z'^-r^ 
orbitals Q]. The mean number of d-electrons per Mn is 4-x, the Hunds rule coupling 
is believed to be very strong relative to the d-d hopping and the spin-orbit coupling 

so the spins of all of the d-electrons on a given site must be parallel. Three of the 
d electrons go into tightly bound core-like dxy, dxz, dyz orbitals forming a core spin Sf 
of magnitude 3/2, to which the outer shell electron (which may hop from site to site) 
is aligned by the Hund's rule coupling. The Hamiltonian containing this physics is 

Hd-cx = - ^t^La^jha - JuY^i ' dJaa^a/jdia/? (1) 

(ij)a,b,a 

Here dl^^ creates an electron in an outer-shell orbital state a = x^ — y^ or 3z^ — r^ 
and spin a, Jh is the Hunds rule coupling connecting the core spin to the outer shell 
electrons, and the interesting limit, which we shall take, is Jh —>■ oo. To study Eq. (1) 
it is convenient to parametrize by polar angles 6i, 0; and to rotate the electrons so 
that the spin quantization axis at, site i is parallel to on site i and then project on to 
the component parallel to S?. The matrix Ri which accomplishes this is Rj = cos^^i/2 

1 -|- isin(6'i/2)sin0icr^ -|- isin(^i/2)cos0icr ^. The electrons may be integrated out, and 
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the partition function Z written as 

VcoseV(j)i expA (2) 



with the action A given by 

A= Tr In [9, - /i + R+drRi - t|^^(RjRj + RjRi)]ii 

(3) 

+iSc /o^ dr 0i(l - cos6'i) 
Here the first term comes from integrating out the electrons and the second is the 
Berry phase term for the core spins. The subscript 11 on the argument of the loga- 
rithm comes from the requirement that in the Jh ^ cxo limit the outer shell electron 
must be parallel to Sc. At low T we may expand about the ordered ferromagnetic 
state 6*1 = 0. The effective action becomes A = Af + Agw Ap is the free energy 
for free fermions moving in the band structure defined by tf^. Asw may be written 
in terms of the magnetization variables Mx and My describing deviations from the 
ordered state (with magnetization taken parallel to z) as 

Aw = I lo dr / i (S, + i^) Mk X 9.M_k 
-2Ka2k2Mk ■ M_k + C(M^) 

Here K = J2ah tfi+i(cJaCi+ib) is the electron stress- energy tensor and is related to the 
integral of the optical conductivity as described below, a is the lattice constant. Here 
c} creates a spin polarized electron on site i. Equation (4) is the action for a quantum 
ferromagnet with spin S* = Sc + and stiffness K. This action implies that the 
magnon dispersion is 

c.^mag = (K/S*)(ka)2 (5) 

A very similar result for the magnon dispersion was obtained by Kubo and Ohata 
using different methods. 

The quantity K is very important because it is the only energy scale in the theory. 
We may estimate K ~ 2tn, where n is the electron density and t the hopping energy. 



A recent band theory calculation found a bandwidth of 2eV implying t ~ 0.25eV and 
K ~ O.leV for n = 0.3. |Q|.From K we may, e.g. estimate the ferromagnetic transition 
temperature as follows: the cubic lattice Heisenberg model with exchange constant J 
has a magnon dispersion uj = 2JS(ka)^ ||^. The known relation between J and Tc 
then implies a Tc ~ 2.9K(S* + 1)/S* ~ .3eV, more than an order of magnitude higher 
than the observed value. This discrepancy, we believe, is evidence that additional 
physics, not included in the double exchange model, is important for Lai.xSrxMnOs. 

A direct measurement of the magnon spectrum would determine K. In the ab- 
sence of this measurement, one may estimate K from the optical conductivity a{uj). 
In a one-band model with only nearest neighbor hopping duja{uj) = vre^K/a [pT| . 
To extract K from conductivity data on a real material one must remove the in- 
terband contributions to a. Ambiguities arise because there is often no clear de- 
marcation between interband and intraband contributions. The optical conductivity 
of Lai_xSrxMn03 has been measured for x = 0.175 and x = 0.3 |[12|. Roughly, 
cr{uj) ^ {500fiQ — cm)~^ ^ 0.3 eV, independent of u and it seems reasonable to 
assume that for u < leV the conductivity is dominated by the conduction band. 
Using a = 4A, one finds K ^ 0.03 eV ~ 400K, much less than the band structure 
estimate. Note that even this value of K implies a magnetic transition temperature, 
much higher than the observed Tc ~ 200K for x = 0.175. 

We now turn to the properties of the model in the regime T ~ Tc. Because 
we have already shown that the low T properties are those of a quantum model 
with the relatively large spin value S ~ 2, it seems reasonable to suppose that near 
Tc we may consider classical spins and so neglect the imaginary time dependence 
of the angular variables in Eq. (3). The problem then becomes that of electrons 
of electrons moving on a lattice with hopping amplitude t^*^ = ty'^[cos(^^i/2)cos(^j/2) 
+cos(0i — 0j)sin(6'i/2)sin(^j/2)]. We further assume that contributions to the partition 
function in which fermions move on closed loops in real space may be neglected. We 
may then rotate the 0i independently and therefore replace this by the familiar double 
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exchange form tj^*^ = ty'^y 2S^''^^ [§!• then replace Eq. (1) by 



Hcff = - XI ^ 



ijab v2 



l + ^(cicjb + h.c.) (6) 



where the Sj are now understood to be classical spins. The free energy function de- 
scribing the spin distribution is to be obtained by integrating out the conduction 
electrons. For any fixed distribution of spins the problem is one of conduction elec- 
trons moving in a lattice with random hopping. To a good approximation, the free 
energy of the conduction electrons depends only on the average hopping so the 
spin energy E({Si}) is given by 



E({Si}) = -T^lnfl + e^^^^-'^)] (7) 
k 



with 



ek = — 2t(coskxa + coskya + cosk^a) 



and 



t = (t?\/l + ^) (9) 

In particular, if the temperature is less than the Fermi temperature of the electrons, 

t^^ I Q. . c. 

E({Si}) = - E + ^fe) (10) 

In other words, the spin energy involves nearest neighbor coupling with scale 
again set by the electron kinetic energy. In the nearest neighbor Heisenberg model 
at Tc, (Si ■ Sj)/S^ ~ 1/3 [0, so an expansion in Si • S/S^ is reasonable, and we may 
conclude that the spin energy is given by the nearest neighbor Heisenberg model, 
with J = K/2^/2. Thus thermal effects do not significantly change our estimates 
of the energy scales, which are much too large to explain the observed transition 
temperature. 



We now turn to the resistivity of the model near Tc. Before presenting the details 
of the calculations, we make some general comments. The electron-spin-fiuctuation 
interaction in Eq. (6) leads to an electron self-energy with real and imaginary parts. 
The real part leads to a contribution to the electron velocity which increases as the 
temperature decreases and expresses the physics that as (Si • Sj) increases, so does the 
electron hopping. The imaginary part leads to scattering, due physically to fluctua- 
tions in Si ■ Sj. If the spins are treated classically, this scattering is mathematically 
identical to conventional impurity scattering, and leads to a resistivity proportional 
to (pf^)~^ where pp is the electron Fermi wavevector and ^ is the mean free path. 
The mean free path in this static spin approximation is a purely geometric property 
determined by the amplitude and spatial correlations of the fluctuations in Si ■ Sj and 
in particular is independent of the electron velocity. Therefore to compute the resis- 
tivity it suffices to calculate the scattering of the electrons off the spin fluctuations, 
neglecting the possibly large velocity renormalization. Further, we shall see that the 
scattering is sufficiently weak that the Born approximation suffices. Finally, we note 
that because the scattering is static the optical conductivity must have essentially 
the Drude form. 

There is no rigorous expression for the dc resistivity. To obtain a reasonable 
approximate expression we use the convenient and at least qualitatively accurate 



"memory function" method [jT4[ in which one defines the memory function M(ci;, T) 



via 



oo 

lUjt I 



M(^,T)= / dte'-*([H,j]„[H,j]o) (11) 
Jo 

where the current operator j is 

j = iEtf (4cjb - 4cia)(l + fracSi ■ S^S^) (12) 

ijab 

and H is given by Eq. (6). The Heisenberg representation is assumed and the sub- 
script on the commutator denotes the time argument of the operators. The memory 
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function rigorously determines the leading term in a high-frequency expansion of the 
conductivity. By assuming that this leading term is the first term of an expansion of 



a{uj,T) = Ke^/a ito + ^'■^''^^ one finds ||T4| that the temperature dependent resis- 



tivity p(T) = e'^M{uj = 0, T). We have evaluated M{lj = 0, T) from Eqs. 6 and 11 to 
leading order in 1/S and kpa. We find 

p(T) = e^ ^ (S(0) ■ S(-5i)S(R) • S(R + 52))/S^B(R) (13) 

R,5i,52 

Here R labels sites on the cubic lattice, and 6i and 62 are any of the vectors x, y, z 
connecting a site to one of the nearest neighbors. B is proportional to the electron 
current-current correlation function weighted by a factor accounting for the ineffec- 
tiveness of small q scattering in degrading the current. In the free electron approxi- 
mation in which the fermions have a dispersion and = t^ = t, 
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B(R) 



32(pFa) 



sin^PF(R + x) sin^pplR — x| 2sin^pFR 



(14) 



(Pf|R + x|)2 (Pf|R-x|)2 (pfR)2 

Here pf is the Fermi wavevector. 

It is interesting to compare Eqs. 11 and 12 to the expression for p(T) given by 



Langer and Fisher who considered the general question of resistive anomalies 
at magnetic critical points. They began from a model in which the carrier-spin 
coupling was Hc_=s J2i Si • (?ci where Si is a local moment and cXc the carrier spin 
density at site i, and obtained a formula rather similar to Eqs. 11, 12 except that 
instead of the four spin correlator they obtained simply (S(0) • S(R)) because in their 
model, local fiuctuations of Si scatter the electrons while in the model defined by 
Eq. 6 local fiuctuations of Si ■ Sj scatter the electrons. Langer and Fisher found that 
there are two sorts of resistive anomalies at a ferromagnetic transition: for T > T^, 
dp/dT ~ (T — Tc)~" where a is the specific heat exponent (which is believed to 
be slightly negative for the three dimensional Heisenberg model) while for T < Tc 
there is an additional contribution to the resistivity proportional to the square of the 
magnetization. These general conclusions apply also to the present model. 
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One important consequence of the difference in models is that the peak at T > Tc 
found for small kp by Langer and Fisher is absent in the present model. The point 
is that the function B has range kp^. In the range where this is greater than the 
magnetic correlation length, the two spin correlator in the Langer-Fisher expression 
for p diverges in the same way as the uniform susceptibility, i.e. as In the present 
case, the quantity in square brackets is the long- wavelength fluctuations in Si ■ Sj, 
which has the weaker divergence ^4-d-2r?^ ^^^q calculations we have performed for 
0.5 < kpa < 1 the divergence is not visible. 

For T < Tc, or in a magnetic field the same formula applies except that one must 
distinguish between longitudinal and transverse fluctuations, and one must add terms 
in which two of the spins in Eq. 11 are replaced by the uniform magnetization. These 
terms lead to a contribution to p proportional to which is absent in the model of 
Langer and Fisher. In particular, in the model of Langer and Fisher, the only effect 
of a nonzero M was to decrease the total amplitude of the spin fluctuations, leading 
to a decrease in p below T^. In the present model, this effect competes with the four 
additional scattering terms of form M^(S(0) • S(R)) which lead to an increase in the 
resistivity. 

We have evaluated Eq (11) using Eq. 12 for B and calculating the spin correlator 
in the spherical model. Results are shown in Fig. 1 for ppa = 1 and two magnetic 
fields: H = and H = O.lTc. Results for ppa = 0.5 are very similar. Note that 
the results are consistent with the predictions of Langer and Fisher: in the spherical 
model a = — 1 so C and dp/dT have derivative discontinuities at Tc and below Tc an 
additional term, proportional to M^, is operative. However, in the spherical model 
approximation used here the sign of this term is positive: the resistivity increases 
below Tc or in a field. 

The resistivity implied by Eq. 6 has been previously calculated by Kubo and 
Ohata 1^, Searle and Wang [g], and Furukawa [^]. Searle and Wang and Furukawa 
used mean field approximations in which all spin correlations are neglected, i.e. (Si • 
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Sj) — (Si) (Sj) = 0. It is evident from the previous discussion that these correlations are 
essential. Furukawa used an "infinite dimensional" approximation in which he found 
that for T > Tc the core spins fiuctuated very rapidly (i.e. on the scale set by t) and 
led to an enormous imaginary part (^ Jh) to the electron self energy. It is difficult 
to reconcile these results with those presented here. In our work Jh drops out of the 
problem and the Si are seen to be well described near Tc by a classical Heisenberg 
model which entails fiuctuation rates of order T or less. Kubo and Ohata obtained 
via a different method an expression very similar to our Eq. 18, but evaluated the 
spin correlation function using an approximation which neglected the fluctuations in 
(Si ■ Sj) which are responsible for the up-turn we flnd in p near Tc. 

Some representative experimental data from ref |jT6| are shown in the inset to 



Fig. 1. Similar experimental results have been obtained by many authors [^^|1^|T8 



Although the qualitative temperature dependence calculated for T > Tc is consistent 
with the data several important discrepancies are evident: the calculated resistivity 
has the wrong magnitude (by several orders of magnitude), a far too weak doping 
dependence, and an incorrect behavior for T < Tc or in a fleld. Some of these 
discrepancies may be due to the inadequacy of the spherical model, but the magnitude 
and doping dependence cannot easily be explained away. The results may be traced 
to the fact, evident already in Eqs. 9 and 10, that in the double exchange model the 
magnetic fluctuations are a weak perturbation on the electron hopping and imply a 
kp^ > 1. 

The discrepancy suggests that some other mechanism, not present in the double 
exchange model, must act to substantially reduce the electron hopping. We suggest 
that this mechanism is a polaron effect due to a very strong electron-phonon coupling 
stemming from a Jahn- Teller splitting of the Mn^^ ion. Other authors, most notably 



Kusterset. al. ||17|| have argued in favor of a magnetic polaron picture. Our calculation 
shows that the standard double exchange Hamiltonian does not contain magnetic 
polaron effects because the effective carrier-spin interaction is too weak to lead to the 
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formation of polarons. On the other hand, the Jahn- Teller coupling is very strong. 
It causes the cubic-tetragonal transition observed at T* ^ 800K in LaMnOa 0, and 
in fact, T* is a dramatic underestimate of the basic Jahn- Teller energy. Using the 
standard Jahn- Teller Hamiltonian and the measured [|T^ oxygen displacements one 
finds that this energy is ~ 1 eV, much greater than the measured electron kinetic 
energy at x = 0.175. It therefore seems likely that the Jahn- Teller energy remains 
important even in the metallic regime 0.2 < x < .45. In this picture the physics 
would involve a crossover between a high T, polaron dominated disordered regime 
and a low T metallic ordered regime. Understanding this crossover requires a theory 
of the interplay of polaron and metallic physics which is beyond the scope of this 
paper. However, it is clear that if polaron physics reduces the mean free path to less 
than a lattice constant, the argument previously given that the carrier resistivity is 
independent of the mass does not apply. Mathematically, if the electron self-energy 
is very large, than the hopping part of the electron Green function is proportional to 
t, not 1/t and so the memory function, Eq. (11) scales as t*-^-*, not t*-'^-'. Then the 
increase in velocity for T < Tc or H 7^ will compete with the extra scattering terms 
and may lead to drop in resistivity for T < Tc, as observed. 

In conclusion, we have presented and compared to data a solution of the "double 
exchange" Hamiltonian widely believed to describe the physics of Lai_xSrxMn03. 
We noted the existence of a relation between an optical property (the low frequency 
spectral weight) and a magnetic property (the spin- wave stiffness). Experimental 
measurements of the spin wave stiffness would be a useful test of the model. We 
showed that the calculated resistivity is much too small, and has an incorrect field 
and temperature dependence, and we proposed that polaron effects are responsible 
for the discrepancy. 
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FIGURES 

FIG. 1. Resistivity calculated from double exchange model as described in the text. 
The solid line is the resistivity in zero field; the dashed line is the resistivity in a field of 
0.15Tc. The inset displays data from Tokura et. al.. 
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